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has been proposed for Yang-Mills theories, inspired by an analogy with disordered systems in 
condensed matter physics. This is based on a particular average over Gribov copies which avoids 
the Neuberger zero problem of the standard Fadeev-Popov construction. The proposed gauge 
fixing can be formulated as a local renormalizable field theory in four dimensions and is well- 
suited for analytical calculations. A remarkable feature is that, for what concerns the calculation 
of ghost and gauge field correlators, the gauged-fixed action is perturbatively equivalent to a 
simple massive extension of the Faddeev-Popov action. The renormalization group flow of the 
theory admits infrared safe trajectories, with no Landau pole. The one-loop calculations of Yang- 
Mills two-point correlators show remarkable agreement with lattice simulations all the way to the 
deep infrared. 
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1. Introduction 



A proper formulation of gauge fixing in non-abelian gauge theories is a longstanding issue. 
The existence of Gribov copies for the most common choices of gauges render the standard Fadeev- 
Popov (FP) construction problematic at least away from the perturbative regime [|0. For instance, 
in these gauges the FP construction for a SU(N) Yang-Mills (YM) theory discretized on a finite 
lattice is plagued by the so-called Neuberger zero problem, inherent to any BRST invariant gauge 
fixing [0]. In the following, we consider the Euclidean theory in the Landau gauge d^A^ = 0, which 
has been the most studied on the lattice. In that case, Gribov copies corresponds to the extrema of 
the functional &[A] = { v tr{A 2 } along the gauge orbit A u of a given field configuration A. 

A way to cope with the Gribov issue is to pick up only one copy for each field configuration. 
For instance, finding a minimum of & [A u ] is a relatively easy task in lattice simulations [Q. How- 
ever, it is not known how to implement this procedure for continuum approaches, which renders 
comparisons somewhat tricky [Bj. Gribov and Zwanziger (GZ) |^[ |6|] have proposed to restrict the 
path integral over gauge fields to the first Gribov region, where the FP operator has only positive 
eigenvalues. This can be approximately formulated in terms of a local, renormalizable field theory 
[|7|] at the price of introducing a collection of auxiliary fields and a dimensionful parameter which 
controls the (soft) breaking of BRST symmetry. However, the first Gribov region is not exempt 
of Gribov copies: ^[A u ] has in fact many minima. A refined version of this proposal, where one 
introduces phenomenological vacuum condensates, agrees well with lattice data 

In |]l]|, we have proposed a different approach to the Gribov problem based on averaging over 
Gribov copies, instead of trying to single out a particular one, in a way that can be formulated 
in terms of a local action. The averaging weight can be chosen so as to break the nilpotent BRST 
symmetry soflty, leading to a perturbatively renormalizable gauged-fixed theory in four dimensions. 
Furthermore, a non-flat weight lifts the degeneracy between copies and avoids the Neuberger zero 
problem. The proposal of [jl]| is very much inspired from an analogy between the Gribov problem 
and that of dealing with potentials presenting a landscape with (exponentially) large number of 
nearly degenerate minima in the physics of disordered systems 

In practice, we use the functional ^[A u ] to weight the various copies. This particular choice 
allows for an elegant and powerful formulation in terms of a collection of replicated supersym- 
metric gauged non-linear sigma (NLa) models. A remarkable consequence of the symmetries of 
the theory is that the NLa model sector of the theory essentially decouples in perturbative calcula- 
tions of YM ghost and gluon correlators. In that case, our gauge-fixing procedure is perturbatively 
equivalent to a simple massive extension of the FP action in the Landau gauge, which is a par- 
ticular limit of the Curci-Ferrari (CF) model [|RJ. This provides a fundamental link between this 
phenomenological model and YM theories. Recent one-loop calculations of two-point correlators 
in the CF model have been shown to agree remarkably well with lattice data all the way to the 
deep infrared [11]. Moreover, the presence of a new running parameter (the gluon mass) opens the 
possibility for infrared safe renormalization schemes, with no Landau pole, as discussed in [ 1 1]. In 
this contribution, I briefly review the main aspects of the proposal of [jl]], recalling some one-loop 
results of [11]. I discuss some open issues and mention ongoing research in this context. 
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2. Weighting Gribov copies 

We consider a SU(N) YM theory on M. d , with classical action 

S YM [A] = ^fti{F 2 }, (2.1) 

where F^ v = d^A v — d v A^ — igo[A^,A v ] and J x = J d d x. For any operator ff[A], we define the 
average over the Gribov copies A u ® of a given field configuration A as: 

(ff\A\) = Ll 1 l{) v , (2.2) 

where s(i) is the sign of the functional determinant of the FP operator —d^D^ evaluated at A = A u< & 
and 5w [A] is a given positive definite weight functional. Here, the sum runs over all Gribov copies 
of the Landau gauge, that is over all extrema of ^[A u ] for a given A. Eq. ( ]22| ) defines a good 
gauge fixing in the sense that it does not affect gauge-invariant operators: (^^[A]) = <^inv[A] for 



[A ] = ^jnvfA]. The denominator in (2.2) is crucial for this property to hold. In the case of 



a flat weight, Sw[A u <- i '>] = const, all copies are degenerate -as in the case of the FP construction- 
and this denominator -as well as the numerator in the case of gauge-invariant operators- vanishes: 
= 0- F° r a non-flat weight the degeneracy is lifted and there is no Neuberger zero problem. 
In [jl|], we choose to weight copies with the functional &\A\ itself: 

S w [A] = Po ftv{A 2 }. (2.3) 

J X 

with j3o > 0. This interpolates between the usual FP construction for j8 — > (see below) and the 



absolute Landau gauge [12], which corresponds to selecting the absolute minimum of JF[A ], in 



the limit j8q — > °°. We also note that, assuming a gap between the minima of &[A ] and the first 



saddles, the weight (2.3) suppresses copies outside the first Gribov region for not too small jSq. 



Once the average ( |2.2P over Gribov copies has been performed for each individual field con- 
figuration one performs the usual average over YM field configurations, hereafter denoted by an 
overall bar. The average of a given operator G\A\ in our gauge-fixing procedure is thus obtained as 
a two-step average 1 : 

f®A(0[A])e- S ™ [A \ (24) 



/^Ae" s YM[A] 
2.1 Field theoretical formulation 



The discrete sums over Gribov copies in ( 2.2 ) can be written as constrained functional integrals 



over local elements U (x) of the gauge group. The constraint d^A 1 ^ = can be exponentiated by 
means of a Nakanishi-Lautrup field h. Similarly, the corresponding Jacobian multiplied by the sign 
s(i) in ( |2~2l ) is nothing but the determinant of the FP operator, which can be exponentiated by means 
of standard FP ghost fields c,c. Denoting collectively "V = (U,c,c,h), we have: 



'A somewhat similar gauge-fixing has been proposed in [^] where, however, the average was not restricted to 
Gribov copies in the Landau gauge. This difference is essential, e.g., in making the present proposal renormalizable. 
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with the gauge fixing action 

S GF [A, t) = ^tr{j3oA 2 + 2d M cD M c + lihd^ }^ . (2.6) 

The latter presents a number of linear and non linear symmetries, including a generalized non- 
nilpotent BRST symmetry. These are most easily seen by introducing the matrix superfield 

6,6) = jgotfetfe+eeQu (27) 

living on a superspace made of the original Euclidean space W l supplemented by two Grass- 
mann dimensions (6,6), which we collectively denote by (9. Here, h = ih — i^{c,c} and the 
^-dependence only appears through the fields U,c, c and h. It is straightforward to show that 

Sgf [A , V } = J ^ tr { {% V ) 1 ( % V ) } , (2. 8) 

where we introduced the covariant derivative = duY + igtfVA^ and f e = J d6d6 g l / 2 (&), 

withg 1 /2(e) = (j3 o 00-l) , can be seen as the invariant measure associated with curved Grassman 
dimensions JT4|]. Eq. ( |2.8[ ) is the action of a supersymmetric gauged NLa model on a curved super- 
space. It is invariant under the isometries of the curved superspace, the super gauge transformation 
A M -)• Af = ^A^- 1 + j-^d^- 1 and f -> r^" 1 , where ^ = W(x,6) is a local element 
of SU(N) on the superspace, as well as under the right SU(N) transformation "V — > ^rV , where 
the matrix = %r(6) can be local in Grassmann variables. 



A non-trivial issue in Eq. fl2.4| ) is the presence of the denominator of the average fl2.5[ ) over 
Gribov copies, which is a highly non-linear, non-local functional of the gauge field A. Similar 
two-step averagings are common in the theory of disordered system and are efficiently dealt with 



by means of the replica trick []15|]. In its simpler version, the latter amounts to writing, formally, 

.n-l 



1 



lim [V\ (^fke-^^A. (2.9) 

n-^OJ ^ V / 



j c$y e -5bF[A,-r] 

Introducing a «-th replica from the numerator in ( |2.4[ ), the final average over gauge fields then reads 

where S^A,!^}] = Sym[A] + Y!l=\ ^gf[A, For perturbative calculations, one needs to factor 
out the volume of the gauge group. This can be done by exploiting the gauge invariance of the 
integration measure %■ A and of the YM action and using appropriate changes of variables to extract, 
say, a factor / S>U n . Renaming (c n ,c n ,h n ) — > (c,c,h), one finally gets 

— ss(A,c,c,h,m)me-^ My}] r2tn 

"n™ j&(A,c,c,h,{r})e- s i A ^Mn] ' K ■ ' 



with &(A,c,c,h,{r}) = &(A,c,c,h) (U" k =\^k) and 

S[A,c,c,h,{r}} = f tr{iF 2 + j 8oA 2 + 2^cD M c + 2/HA4 + AE / tr \{%^kf {%%)}■ 



(2.12) 
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This describes a collection of n — 1 gauged supersymmetric NLa models coupled to a gauge-fixed 
YM field with gauge fixing Sgf[A,1,c,c,/i]. Notice that, as is clear from Eqs. ( |2.6| ) and (2.8). 
each replica produces a term /3oA 2 , thus giving rise to a bare mass term raj3oA 2 . Thanks to its large 
number of symmetries, the theoiy ( 2.12 ) has been shown to be perturbatively renrormalizable in 
d = 4, with only two renormalization factors, in [fil]. 



3. Perturbative equivalence with the Landau limit of the Curci-Ferrari model 



Perturbation theory is most conveniently formulated in the supersymmetric formalism, which 
makes transparent the (dramatic) consequences of the supersymmetries for loop diagrams. Para- 
metrizing the constrained superfields % in terms of unconstrained ones, e.g., % = exp{/goA*}, it 



is straightforward to obtain the various propagators of the theory ( |2.12| ), written below in Euclidean 
momentum space. The gluon propagator reads: 



A a ,{ P )A h v (-p) 



lab 



p 2 + np 



PfiPv 



(3.1) 



where the square brackets represent averaging with the action ( |2.12[ ) with n ^ 0. It is transverse in 
momentum, as a result of Landau gauge, and massive, with bare square mass nfio, as a result of our 
particular gauge fixing procedure. The ghost propagator assumes the usual form: 



c a { P )c h (-p)\ =8 ab /p 2 . 
Finally, the superfield propagator is given by 

A2(p,e)A?(- jP> e , )l = s ab 8 kI d(e,e!)/p 2 ., 



(3.2) 



(3.3) 



where 8(6_,6^) = g~ l ^ 2 (d_) (6 — 6 )(8 — 6') is the Dirac function on the curved Grassman space. 

The vertices of the action ( 2.12j ) which do not involve the superfields A* are the same as for 
the usual FP Landau gauge. The NLa model sector of the theory gives vertices with an arbitrary 
number of legs and either zero or one gluon leg. Clearly, the latter are local in Grassmann 



variables. Using ( |3.3| ), we conclude that any closed loop of A^ superfields with p vertices insertions 
involves the product 8(6_i,d_ 2 ) " " " S{6_ pl d_i) °= S^,^) = 0. Hence, the NLa model sector of the 
theory is tree-level exact. This is not surprising since, in the above construction, the role of these 



superfields is in fact to reconstruct the weighted sum (2.2) on Gribov copies, that is to project on 



the extrema of the action (2.3) 



This observation has two important consequences. First, the only perturbative source of de- 
pendence in the number n of replicas is the bare gluon mass in ([3j]). Second, correlators or vertex 
functions involving only the fields A, c and c do not receive any loop contribution from superfields. 
They are thus obtained from the very same diagrams as in the FP Landau gauge, with usual YM 



vertices and with propagators given by (gjj) and ( |3.2| ), that is from the effective action 



S\A,c,c,h,{n) 



S en [A,c,c,h} = Jtvi^F 2 + n^)A 2 + 2d tl cD ll c + 2ihd ll A ll ^ . (3.4) 



This simple massive extension of the FP action is a particular case of the CF model [1C]. 
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A key point here concerns the issue of the limit n — > versus renormalization. It is clear from 
the above considerations that taking first the zero replica limit leads to the standard FP Landau 
gauge. This, however, is not satisfactory since one expects gauge-dependent quantities to depend 
on the gauge-fixing parameter j3o- A similar situation arises if one introduces a renormalized gauge- 
fixing parameter, e.g., as j8o = Zp(3. An alternative renormalization scheme is to redefine the square 
mass as nj8o = Z m 2tn 2 . This absorbs the remaining ^-dependence and gives a non-trivial n — > limit. 



4. One-loop results 

The ghost and gluon propagators have been computed at one-loop in the theory ( |3.4[ ) in [p"T|]. 
The results are in remarkably good quantitative agreement with lattice data for SU(2) and SU(3) in 
d = 4 and give a fairly good qualitative description in d = 3. Fig. [j] reproduces the gluon and ghost 
propagators Gf v (p) = 8 ab (8^ v - PilPv / p 2 )G{p), and D ah {p) = 8 ab D(p) computed in [0]. 




p (GeV) p (GeV) 

Figure 1: The gluon and ghost dressing functions, p 2 G(p) and F(p) = p 2 D(p), for the SU(3) theory in 
d = 4, obtained in [11]. Lattice results are shown as blue dots. The lines show the one-loop results obtained 
from the action (3.4), without (red) and with (black dashes) renormalization group improvement in the 
UV. The employed renormalization conditions are G~'(0) = m 2 , G _1 (/i) = m 2 + p 2 , D~ l (p) = p 2 and a 
standard Taylor scheme for the coupling. The parameters are m = 0.54 GeV and g = 4.9 at p = 1 GeV. 

Another interesting observation is that the theory ( pM| ) admits infrared safe renormalization 
schemes, with no Landau pole, as discussed in [ [TTj , p. 2 Defining the renormalized fields and 
constants A = \fZ\A r , c = \fZ c c r , c = \/Z^c r , go = Z g g and n% = Z m 2tn 2 , the following renormal- 
ization conditions 



G- l (p)=m 2 + p 2 , D-\p) = p 2 , Z A Z c Z m 2 = \, (4.1) 

together with the Taylor condition for the coupling, Z g Z c \fZ~A = 1, lead to the renormalization 
group (RG) flow depicted in Fig. |2[ There is an ultraviolet (UV) attractive fixed point at m = 
and g = 0: both the bare coupling go and the bare mass riffa vanish in the process of removing 
the UV regulator. We note that this is compatible with taking the limit n — > at fixed gauge-fixing 
parameter j3o- Most remarkably, RG trajectories fall in two distinct classes, depending on the initial 
conditions in the UV: those with or without a Landau pole in the IR. It is worth mentioning that the 
best values of m and g for describing lattice data in the scheme (LI ) belong to the second class [pT|], 



2 The IR flow of this theory has also been studied in Jl6| ] 
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Figure 2: One-loop RG flow in the plane (in — m/fi,g). The arrows indicate the flow towards the infrared. 



5. Issues and perspectives 

The gauge-fixing proposed in ^ provides what appears to be an essential feature of the gluon 
propagator, an effective mass term, already at tree level. This raises the following issues. In prin- 
ciple, different values of the gauge-fixing parameter j3o, corresponding to different RG trajectories 
in Fig. ||, could give (vastly) different results for the YM correlators. Instead, lattice results in the 
minimal Landau gauge show at best a mild sensitivity with the selected Gribov copy [|7|. However, 
it is conceivable, as mentioned previously, that the minima of ^[A u ] be nearly degenerate and well 
separated from the first saddles. This is supported by numerical investigations on small lattices, 
where all copies can actually be found [|T|]. If so, there exists some range of values of j8o which 
corresponds to giving essentially an equal weight to copies in the first Gribov region -probed by 
lattice simulations- and to suppressing those outside the first region. It is presently not known how 
lattice results are affected if one selects a copy outside the first Gribov region. This issue requires 



detailed studies of Gribov copies, including sadles, in the spirit of Ql8p. 

Another important question concerns the limit n — > versus renormalization which, as already 
mentioned, do not commute. A similar issue arises in the theory of disordered systems, where 
one is concerned with the thermodynamic limit instead of renormalization. In this context, it is 



understood that the limit n — > should always be taken last [|15|]. In fact, replicas can be viewed as 
non-trivial external fields which allow one to probe the complicated landscape of the potential and 
which should be removed only at the end of the calculations in order to probe non trivial physics, 
such as spontaneous symmetry breaking. We believe this picture can be fruitful in the YM context 
too. Sending n — > naively leads to the standard FP Landau gauge with nilpotent BRST symmetry. 
Taking, instead, the limit n — > after having properly removed the UV regulator, one captures non- 
trivial IR physics where, however, the nilpotent BRST symmetry is broken. This is very transparent 



after the superfields have been integrated out in (3.4). The mass term «j8o can be seen as an external 



source which explicitly breaks the BRST symmetry. The fact that the symmetry is not recovered 
after this source is eventually removed signals that it is spontaneously broken by IR fluctuations. 3 
To conclude, we believe the proposal of [[Tj] opens a possibility to access non-trivial IR physics 
by perturbative methods. Ongoing research in this context includes the calculation of higher-order 
corrections and higher-order vertices in the theory ( $A^ and the inclusion of quarks [|l9|], or the 



calculation of YM correlators and thermodynamics at finite temperature [20]. It may be of interest 



3 We mention that a similar construction can be made for QED. In this case, however, the mass term does receive 
nontrivial renormalization and thus vanishes in the limit n — > 0. 
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to investigate the non-perturbative aspects of the theory ( [2. 12 ) either with continuum methods or 



with lattice simulations, in the spirit of pl[], see also Other interesting questions concern the 
the generalization of the approach of [|]] to other gauges $3S\, or the relation with other proposals 



such as, e.g., [|23|]. Finally, a major open question is that of unitarity [g4j]. 
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